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Abstract: The following problem, which stems from the "flux phase" problem in condensed 
matter physics, is analyzed and extended here: One is given a planar graph (or lattice) with 
prescribed vertices, edges and a weight \t xy \ on each edge (x,y). The flux phase problem (which 
we partially solve) is to find the real phase function on the edges, 9(x,y), so that the matrix 
T := {|t a; j / |exp[^(x, y)]} minimizes the sum of the negative eigenvalues of —T. One extension 
of this problem which is also partially solved is the analogous question for the Falicov-Kimball 
model. There one replaces the matrix —T by —T + V, where V is a diagonal matrix representing 
a potential. Another extension of this problem, which we solve completely for planar, bipartite 
graphs, is to maximize |det T\. Our analysis of this determinant problem is closely connected with 
Kasteleyn's 1961 theorem (for arbitrary planar graphs) and, indeed, yields an alternate, and we 
believe more transparent proof of it. 
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I. INTRODUCTION 

The genesis of this paper was an attempt to understand a problem in condensed matter 
physics related to questions about electron correlations, superconductivity and electron-magnetic 
field interactions. The basic idea, which was proposed a few years ago is that a magnetic field can 
lower the energy of electrons when the electron density is not small. Certain very specific and very 
interesting mathematical conjectures about eigenvalues of the Laplacian were made and the present 
paper contains a proof of some of them. Furthermore, those conjectures lead to additional natural 
conjectures about determinants of Laplacians which we both present and prove here. It is not clear 
whether these determinantal theorems have physical applications but they might, conceivably in 
the context of quantum field theory. Some, but not all, of the results given here were announced 
earlier in [LE]. 

The setting is quantum mechanics on a graph or lattice. (All our terminology will be defined 
precisely in the sequel.) Physically, the vertices of our graph A can be thought of either as a 
discretization of space (i.e. replace the Laplacian by a finite difference operator) or they can be 
seen as locations of atoms in a solid. There are |A| vertices. In the atomic interpretation the edges 
become electron bonds joining the atoms, and the model is known as the tight-binding model or 
Huckel model. The natural Laplacian C associated with A is a |A| x |A| matrix indexed by the 
vertices of A and whose diagonal elements satisfy — C xx = number of attached edges (or valency) of 
vertex x. The other elements are C xy = 1 if x and y are connected by an edge, and zero otherwise. 

For us it is more convenient to consider the matrix C which is the Laplacian without the 
diagonal term, i.e., C xx is replaced by zero. In the context of graph theory C is also known as the 
adjacency matrix. There are three excuses for this: (i) in the solid state context C is the natural 
object because atoms do not bond to themselves; (ii) most of the graphs that are considered in the 
physics literature have constant valency, so C and C have the same spectrum modulo a constant 
which is equal to this valency; (hi) mathematically C seems to be the more natural object — from 
our point of view, at least — because its spectrum on a bipartite graph is always a union of pairs 
A and —A (when A 7^ 0), as explained in Sect. II. The spectrum of C generally does not have any 
such symmetry. 

We label the eigenvalues of C by Ai(£) > A2(£) > The Hamiltonian for a single electron 

is —C or —£, and we take it to be —C here. If our system has M free electrons the rule of quantum 
mechanics is that the eigenvalues of our system are all the possible sums of M of the (— Aj)'s in 
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which each — Aj is allowed to appear at most twice in the sum. There are ( 2 j^') eigenvalues. In 
particular, if M = 2N the smallest eigenvalue is 

Ei N) =-2j2^)- (1-1) 

A (spatially varying) magnetic field is now added to the system in the following way. C xy is 
replaced by T xy = C xy ex.p[i6(x, y)], with 9 real and with 9(x, y) = —0(y, x) so that T is Hermitian. 
The function 6(x, y) is interpreted physically as the integral of a magnetic vector potential from 
the point x to the point y. This T is the discrete analogue of replacing the Laplacian on IR™ by 
(V — iA{x)) 2 (with V = gradient), which is the Laplacian on a U(\) bundle. 

The central question that we address is this: What choice of 9 minimizes Eq N ^ for a given N? 

In order to appreciate this question, consider the N = 1 case. Then, 9 = is an answer 
because (with <f> being the normalized largest eigenvector of T) \i(T) = ^ <f> x <f> y L xy exp[i9(x, y)] 
< l&xW&yl^xy < Ai(jC). This proof that 9 = is optimum also works in a more general setting, 
namely for the lowest eigenvalue of the "Schrodinger operator" —T + V, where V is any real 
diagonal matrix. Again, T = £, or 9 = 0, minimizes — \\(T — V). The same is true in IR n for 
— (V — iA{x)) 2 + V(x); the minimum occurs when A(x) = 0. This conclusion is known as the 
diamagnetic inequality and states, physically, that "a magnetic field raises the energy" . 

It was discovered by [AM] and [KG] that the situation can be quite different when is close 
to ||A|. (When N = |A|,£'J Ar) = TrT = for all 9; hence N = ||A| is the most extreme case.) 
Since then, the problem has been investigated for various lattices and A's by several authors such 
as [BR], [BBR], [HLRW], [RD], [WP], [WWZ] some of whom consider it to be important in the 
theory of high temperature superconductivity. [HLRW], for example, start with the square lattice 
Z 2 , take A to be a large rectangular subset of Z 2 , and then let |A| — > oo and N — > oo with A/|A| 
fixed. They also take the magnetic flux (which is the sum of the 0's around the edges of a face, 
and which is defined in Sect. II) to have the same value in each square box of Z 2 . On the basis of 
their numerical evidence they proposed that flux/box = 27rA/|A| is the optimal choice. In [AM] 
the term "flux phase" was introduced to describe this state in which the presence of a magnetic 
field lowers the energy. 

It should be pointed out that the spectrum of £ for Z 2 as a function of constant flux/box was 
discussed by many authors for many years; it was Hofstadter [HD] who grasped the full beauty 
of this object — which is anything but a continuous function of the flux and which is full of 
gaps — and called it a "butterfly". The spectrum can be found by solving a one-dimensional 
difference equation, due to Harper [HP], which is a discrete analogue of , but more complicated 
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than, Mathieu's equation. The spectrum is such a complicated function of the flux that it is 
difficult to decide on the optimum flux for a given N. 

The most striking case is JV/|A| = \, or M = |A|, which is called the half-filled band. The 
optimal flux is supposed to be ir, which is the maximum possible flux since flux is determined only 
modulo 2tt and since flux and —flux yield identical spectra. It is this case that we investigate in 
this paper in an attempt to verify the rule just stated and which appears in [AM], [HLRW], [RD]. 
We are completely successful only in some special cases, but we have been able to generalize the 
problem in several interesting directions. For example, one of our main results is Theorem 3.1. It 
completely solves the problem for determinants (i.e. for products of eigenvalues instead of sums of 
eigenvalues) on bipartite planar graphs. 

Our determinant theorem turns out to be closely related to Kasteleyn's famous 1961 theorem 
about planar graphs, which allowed him to solve (in principle) the dimer problem and Ising model 
for all planar graphs. Our route, via fluxes, gives an alternative proof of Kasteleyn's theorem and, 
we believe, a more transparent one. This is presented in the Appendix. 

The setting we adopt is a general graph A, with no particular symmetry such as Z 2 enjoys, 
and an arbitrary, but fixed amplitude \t xy \ > given on each edge (\t xy \ = 1 in the case of C ). 
The problem is to determine and T := {t xy } x . y& \ with t xy = \t xy \ exp[i6(x,y)} so as to minimize 
the (absolute) ground state energy 

E (T) := — Tr|T| - TrT = — Tr|T|, (1.2) 

with Tr = Trace. The right side of (1.2) is twice the sum of the negative eigenvalues of — T. For a 
bipartite graph this is the sum of the |A|/2 or (|A| — l)/2 lowest eigenvalues of —T. 

A word has to be said here about different definitions of ground state energy. Electrons have 
two spin states available to them and the Pauli exclusion principle states that each eigenstate of 
—T can be occupied by at most one electron of each kind. Thus, each eigenstate can be occupied by 
or 1 (twice) or 2 electrons. That explains the factor of 2 in (1.1): there the lowest N eigenstates 
of —T are each occupied by two electrons. Our definition of Eq(T) in (1.2) is the absolutely lowest 
ground state energy and corresponds to the electron number being twice the number of negative 
eigenvalues. On the other hand, the half-filled band would have the electron number equal to | A| by 
definition. If |A| = 2N is even, then the half-filled band ground state energy is given by (1.1) with 
N = |A|/2. If |A| =2N + 1, the half-filled band ground state energy is -2 ^f =1 Xj - Xn+i- It is 
this half-filled band energy that is mostly considered in the physics literature. However, we regard 
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our definition (1.2) as mathematically more natural and physically as interesting as the strict half- 
filled band definition. For bipartite graphs E (T) and with N = |A|/2 or (|A| — l)/2 agree 
with each other. (Note that if A is bipartite and |A| = 2A" + 1 then Aat+i = 0.) 

The two definitions can produce strikingly different conclusions, however, in special cases. In 
[RD] the ground state energy E (T) of A" electrons (including spin) hopping on a ring of A" sites 
is considered. By Theorem 4.1 we know that for a ring with A" odd (and which is therefore not 
bipartite), the expression (1.2) is minimized by flux tt and flux 0. However it has been shown in 
some cases (see [RD]) that the half-filled band energy for such a ring is minimized by the flux tt/2 
(which, incidentally, we call the canonical flux in this paper). 

There is an important difference between our minimization problem and the one in [HLRW] 
and some other papers in the physics literature. For a regular structure like Z 2 we allow different 
fluxes in different boxes. In the physics literature the problem is sometimes stated with constant 
fluxes or with periodic fluxes. We find our formulation (with arbitrary fluxes) to be more natural 
mathematically and we believe it to be more natural in those physical problems where this theory 
might be applicable. 

Besides the ground state energy problem we consider other functions of T, such as In | det T\ = 
Trln|T|. A particularly important one, physically, is InS where S is the grand canonical par- 
tition function with chemical potential /i and inverse temperature j3, given by 



E = £-"EE"'E ex p 

mi m| A | m n jA 



|A| 

& X] A? ( n i +m,j) + 0n 



|A| 

n{l + exp[/3(A J+ ^)]} 2 , (1.3) 
j=i 



where the sum on each rij and m; is over the set {0, 1}. The physical free energy is defined by 
T = —f3~ x lnS. We consider only jjl = here because that corresponds to a half-filled band in the 
bipartite case (see (4.5) and footnote). 

Another important quantity is the gap, G(T), which is not defined by a trace. We define it 
to be 

G(T) = -X N+1 + X N , (1.4) 

where A" is the number of negative eigenalues of — T, i.e., the smallest number such that Eq N ^ = 
E (T). Clearly G(T) is the energy needed to add one more electron to the system from the 
absolute ground state. For the half-filled band on a bipartite graph with |A| = 2N, G(T) = 2Ajv- 
This, however, may not be mathematically interesting because \n may be automatically zero for 
dimensional reasons. That is, if \A\ and \B\ are the two subsets of vertices of A that define the 
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bipartite structure, then T always has at least \\B\ — \ A\\ zero eigenvalues. For this reason we 
define G(T) for a bipartite A (with |A| odd or even) to be 

G{T) = A| A , - A|B| +1 , (1.5) 

assuming \B\ > \A\. We can then ask the question: Which flux maximizes G(T) or G(T) in the 
bipartite case? 

So far we have discussed free — or noninteracting — electrons. The same questions can be 
asked for interacting electrons and very much less is known in that case. In Sect. VIII, however, 
we are able to carry over our techniques to one example — the Falicov-Kimball model. 

Many of these results were announced in [LE]. We thank P. Wiegmann for bringing this 
problem to our attention and along with I. Affleck, D. Arovas, J. Bellissard and J. Conway, for 
helpful discussions. 

II. DEFINITIONS AND PROPERTIES OF FLUXES 

A graph A is a finite set of vertices (or sites), usually denoted by lower case roman letters 
x,y,z etc., together with edges (or bonds), which are certain unordered pairs of distinct sites 
and are denoted by (x,y), equivalently (y,x). Thus there will be at most one edge between two 
vertices. The set of sites or vertices will be denoted by V or V(A) and the number of them by |A|. 
The set of edges will be denoted by E or E(A). If (x, y) G E the sites x and y are said to be end 
points of the edge (x,y). 

A graph A is connected if for every pair of sites x and y there is a path P in A connecting x 
and y, i.e., there is a sequence of points x = xo,x\,X2, ■■■ ,x n = y such that (xi,Xi+i) is an edge 
for every < % < n. Although A is not just the set of vertices, but also contains the edges, we 
shall nevertheless sometimes write x € A where x is a site in A. 

A hopping matrix T associated with a graph A is a Hermitian |A| x |A| matrix indexed by 
the sites of A, with elements denoted by t xy = t yx for x, y 6 A, and with the important property 
that t xy ^ only if (x,y) € E, i.e., if x and y are connected by an edge. In particular, t xx = 
for all x € V. The T matrix is the important object here. For that reason if t xy = for any 
edge (x, y) we might as well delete this edge from the graph A. Thus, without loss of generality 
we can assume that every t xy is nonzero and that the corresponding graph A is connected. If it 
is not connected T breaks up into blocks which can be considered separately. We call \t xy \ the 
hopping amplitudes. No other assumption is made about t xy unless explicitly stated otherwise. 
The eigenvalues of T are usually denoted by A, and sometimes by A(T) to be more specific. 
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A circuit C of length £ in A is an ordered sequence of distinct sites xi,...,xg with the property 
that (xi,Xi + i) is an edge for % = 1, . . . ,£ with x^+i = X\. We explicitly include 1 = 2. Note that 
X2, ■ ■ ■ , xe, X\ is the same circuit as C, but xg, xg-i, . . . ,X\ is different. 

If C is a circuit then we can define the flux <&c of T through C, which is a number in [0, 2tt), 

as 



$c = arg JJ )= ar § II T ■ 



\i=l / \ C / 

The symbol H T has an evident meaning. 

c 

A gauge transformation is a diagonal unitary transformation U with elements, -u^ = 
ex.p[i(f) x ]5 X y, where 4> x : V — > IR is a function on the sites. Obviously a gauge transformation 
T — > U*TU leaves the spectrum of T and all the fluxes unchanged. 

2.1. LEMMA (Fluxes determine the spectrum). Let T andT' be two hopping matrices 
which have the same hopping amplitudes and the same flux through each circuit C of the graph A. 
Then there is a gauge transformation U such that T' = U*TU . 

Proof: By our convention the t and t' are never zero. Thus W xy = t x /t' satisfies 
iW^yl = 1 for all edges (x,y), and the flux of W through each circuit of A is zero. Let xq be an 
arbitrary but henceforth fixed site in A. For any x we can pick a path P connecting xq and x and 
define (fi x = arg I FJ W ] . The value of 4> x does not depend on the choice of the path because, if 



P' is another path connecting xo and x, we have that arg ^fT W j = arg ^fT W J since the flux 
through the circuit given by connecting x to x along P and then connecting x to x along P' is 
zero. If x and y are arbitrary sites in A with (x, y) 6 E, and if we take P x to be a path connecting 
xo to x and P y a path connecting x to y we observe that 1 = n W W xy Y[ W since P x followed by 

(x,y) is a path from xo to y. But this equals ex.p[i(4> x — 4> y )\W xy and hence W xy = exp[— i{<j) x — 4> y )}- 
Thus t' xy = exp[i(fr x ]t xy exp[— i(j) y ], which proves the lemma. I 

Up to now a graph has been regarded as an abstract object consisting of vertices and edges. 
Now we wish to regard graphs as embedded either in HI 2 or K 3 . This means that the sites of A 
can be regarded as distinct fixed points in IR 3 and each edge (x, y) will be identified with exactly 
one piecewise linear curve between x and y € M 3 . It is convenient to exclude the end points x and 
y in the definition of an edge. We require that any one edge does not intersect the other edges or 
the sites. Circuits are then identified with simple, oriented closed curves. 

Obviously any graph can be embedded in IR 3 but only some graphs, called planar graphs, 
can be embedded in IR 2 . It is these graphs that will mostly concern us in this paper. 
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The set of edges E and sites V of a planar graph, regarded as a set in IR 2 , is closed. Its 
complement is therefore open and this complement has a finite number of connected components 
which we label F ,Fi,F 2 , — We define F to be the open set that contains the point at infinity, 
(i.e. the exterior of the graph). The others we call the faces of the graph A. Each face has a 
boundary and this boundary is composed of a subset of the edges and sites of A. For later purposes 
we call elementary circuits those circuits which are entirely contained in the boundary of a single 
face. 

If the graph is planar a circuit, C, of length greater than 2 will have an inside and an outside. 
The interior, which is an open set, is then the union of a certain number of faces, edges and 
vertices called interior faces, interior edges and interior vertices. We denote their numbers 
by /, e, v. We can speak of the orientation of C as being either positive (anticlockwise) or negative 
(clockwise) according as the winding number with respect to a point in its interior is either +1 or 
-1. 

In general, an arbitrary specification of fluxes through the circuits of A may be inconsistent in 
the sense that there may not exist a choice of T with the prescribed fluxes. Some kind of divergence 
- or closedness condition is needed. In two dimension, however, the following lemma shows how 
fluxes can be specified in a consistent way. 

2.2. LEMMA (Construction of phases from fluxes). Let A. be a planar graph and let 
Fi, F 2 , . . . Ff be its faces. Let <I> 2 , . . . , & be any given numbers in [0, 2n). (We call $ J the flux 
through Fj.) Then there is a function 0(x,y) : E(A) — > [0, 2tt) so that the fluxes defined by (2.1) 
with t xy := exp[i9(x,y)\ satisfy $ c = ^interior faces of c& f or ever y positively oriented circuit 
C on A. 

Proof: Pick a point Zj = (zj,Zj) £ IR 2 in each face Fj and consider the vector (really, a 
one-form) field on IR 2 , A(x) = A(x 1 ,x 2 ), with singularities at the Zj's given by 

A(x)= Yl &${x-Zj) 

all faces Fj of A 

with a(x) = (27r)" 1 [(x 1 ) 2 + (x 2 ) 2 ] _1 (-x 2 , x 1 ). For each pair x,y € V(A) with (x,y) £ E(A) define 
y -, 

9(x,y) = J A, i.e., the integral from x to y along the curve representing the edge (x,y). The 

X 

flux $c through any circuit C on A is given by § c A. By Cauchy's integral formula (or Stokes' 
theorem), this integral equals $>c given above. ■ 

Remark: A more practical and direct way to construct phases 9(x, y) satisfying the conclusion 
of Lemma 2.2 is to concentrate the vector field a(x — Zj) along a line. More precisely, let Lj denote 
some semi-infinite line starting from zj and extending to infinity, but which does not intersect any 
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of the sites of A and whose intersection with an edge is always transverse, i.e., nontangent. We 
orient Lj from Zj outwards. For each ordered pair of sites (x,y), with (x,y) £ E, first orient the 
curve representing (x, y) in the direction from x G IR 2 to y 6 IR 2 . Call this oriented curve e. Then 
let [ij(x, y) = S^ ne (±l), where the sum is over all the points of intersection of Lj with this curve 
e and where the + sign (resp. — sign) is taken if the (counterclockwise) angle from Lj to e is less 
than (resp. more than) ir. Finally, we set 9(x,y) = Y,j =1 & fj,j(x,y). 

It is a fact that every planar graph, A, can be triangulated, i.e., that there is a planar graph 
A' with precisely the same vertices as A and whose edges E' (which are now sets of curves) contain 
E, the edge-set of A, and with the property that A' is triangular. The concept of a planar graph 
A' being triangular means that every one of the faces of A' has as its boundary the union of three 
edges and three vertices. We say that A' is a triangulation of A. It is easy to check that the 
three edges must always form a circuit. Note that, in general, a graph A can be triangulated in 
several ways. 

2.3. LEMMA (Number of triangles in a circuit). Let A be a triangular planar graph 
and let C be a circuit in A of length £ > 2. Let f and v denote the number of interior faces and 
interior vertices of C. Then 

£-f + 2v = 2. (2.2) 

For £ = 2, (2.2) is clearly true (with / = v = 0). Therefore we need consider only £ > 2. 

First proof: Let (3 be one of the £ edges in C. This edge must be part of the boundary of 
exactly one of the inner triangles, which we call r. The boundary of r contains 3 edges, f}\, 02, 03- 
There are two cases. 

(a) 0\ = and (32,03 are interior to C. 

(b) 0\ = (3, 02 is an edge of C, (3% is interior to C. 

In case (a) we consider the circuit C whose edges are the same as those of C except that (3 is 
replaced by the two edges (32,03- 111 case (b) we remove 0\ and 02 from C and replace them with 
/3 3 . It is easy to check that £' — f + 2v' = I — f + 2v and that /' = / — 1. By successively removing 
triangles in this way we eventually have only one triangle left, in which case I = 3, / = 1, v = 0. 

Second proof: Euler's formula says that (total number of vertices) + (total number of faces) 
— (total number of edges) = 1. Since £ also equals the number of vertices in C, we have that 
1 = (v + £) + (/) — (£ + e) , where e is the number of interior edges. Each edge in C lies in the 
boundary of precisely one interior triangle, while each interior edge lies in the boundary of two 
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such triangles. Since each triangle has three edges in its boundary, and since C has I edges, we 
have 3/ = 2e + L Therefore e = 3//2 - 1/2 and 1 = v + / - e = v - f/2 + 1/2. ■ 

2.4. COROLLARY (/ is independent of triangulation). Let A 6e an arbitrary planar 
graph and let A 1 denote a triangulation of A. For each circuit C in A i/ie number of triangular 
faces of A' that are interior to C is independent of the triangulation A'. 

Proof: The result follows from (2.2) since i and v do not depend on the chosen triangulation 
A'. ■ 

With the aid of triangulation we can describe the canonical flux distribution for any planar 
graph A. Choose any triangulation A' and place flux n/2 in every triangular face. By Lemma 2.2, 
this defines phases 9(x,y) on E(A') and hence on E(A). A-priori, these phases might depend on 
the triangulation but, by Corollary 2.4, all triangulation give rise to the same set of fluxes through 
the circuits of A. By Lemma 2.1 the 9(x,y) , s are uniquely defined up to a gauge transformation, 
i.e., 6(x,y) — > 9(x,y) + <f> x — 4> y with the function <f> x being the only quantity that might depend 
on the triangulation. Since the flux distribution is invariant under gauge transformations, the 
canonical flux distribution is well defined! 

Of special interest to us are bipartite planar graphs. In general a bipartite graph is a graph 
A whose vertex set V is the union of two disjoint sets A and B with the property that (x, y) is 
never an edge of A if x € A and y £ A or x £ B and y € B. We shall assume \B\ > \A\. If A 
is a planar bipartite graph the canonical flux will always be tt through every elementary square, 
zero through every elementary hexagon etc. However one has to be cautious about this because 
one could have, for instance, a square with vertices a, b, c, d and a fifth vertex g inside the square 
connected by an edge only to a. In this case our rule says that the canonical flux through the 
circuit a, b, c, d is zero and not tt. 

A special feature of bipartite graphs, planar or otherwise, is that the nonzero eigenvalues of 
any hopping matrix T come in opposite pairs, i.e., if A is an eigenvalue of T then so is —A. This 
follows from T = —U*TU where U is the diagonal unitary matrix with +1 on the A-sites and —1 



elements between A and B sites. 

III. DETERMINANTS FOR PLANAR GRAPHS 

One of the main theorems of this paper is Theorem 3.1 about determinants of bipartite graphs, 
and one of the concepts needed there is that of the dimer partition function D(T) of the graph 
A with hopping matrix T. A dimer covering (or matching) of A is a subset {ei, e2, . . . , e n } of E 
such that every site in A is an end point of precisely one of the e^'s. In general, A has many dimer 
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coverings, but it may have none at all. In particular, if |A| is odd or if A is bipartite and \A\ ^ \B\ 
then there are no dimer coverings. 

We define the dimer partition function to be 

dimer coverings i 

where the product is over all the edges ej = (xi,y{) that constitute a particular dimer covering. 
If \t X y\ = 1 then D(T) is just the number of dimer coverings of A. Note that D(T) depends only 
on the |t X j/|'s and is therefore independent of the fluxes. In particular, D(T) is determined by the 
upper triangular array {t xy } x < y . (See the appendix.) 

3.1. THEOREM (Canonical flux counts dimers and maximizes bipartite graph 
determinants). Let A be a planar graph and let \t xy \ be given positive numbers for all edges (x, y) 
in A. For the canonical flux distribution 

detT= (-1)I A I/ 2 D(T) 2 . (3.2) 

If, in addition, A is bipartite the canonical flux distribution maximizes | det T\ among all flux 
distributions. 

Before proving the theorem we make a series of remarks: 

(i) . Unless \A\ = \B\ in the bipartite case, D(T) = and detT = for every choice of flux. 
In the general case, D(T) = unless |A| is even. More generally, we could consider non-bipartite 
graphs with T of the form Tk = (^^* K selfadjoint. This means that edges are 
added between i?-vertices but not between ^4-vertices. It is then an easy exercise in linear algebra 
to prove that detTx = unless > \A\, and that if \B\ = \A\ then det IV is independent of 
K, i.e., detTR- = detT . As an example, start with the simple square, i.e., |A| = 4 and (1,2), 
(2,3), (3,4), (4,1) are the edges. Theorem 3.1 says that the determinant is maximized by flux = tt 
through the square. Now add a diagonal edge (1,3) with some hopping amplitude |*x,3 1 on this new 
edge. We now have a graph that consists of two triangles. The observation just made says that the 
determinant is independent of the individual fluxes through the two triangles and depends only on 
their sum. The canonical flux distribution, which is 7r/2 in each triangle, is optimal, but is by no 
means the unique optimizer. 

(ii) . In the bipartite case, the sign of the determinant as given in (3.2) is correct for any T, 
not just the canonical T. This follows from the A, — A pairing of the eigenvalues which holds for a 
bipartite lattice. 
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We can go further and define the elementary symmetric functions 

e k (T)= U X ^ T ) (3-3) 

ii<i2<---<ik j = l 

for 1 < k < |A| and with eo(T) = 1. They are the coefficients in the characteristic polynomial 

|A| 

det(T + z) = Y j e k (T)z\ A \- k . (3.4) 

fc=0 

Here, the Ai(T)'s are the eigenvalues of T. For the same reason (the A, —A pairing) we see that for 
a bipartite graph 

e fc (T) = 0, k odd 

(3.5) 

e k {T) = {-l) k ' 2 \e k {T)l k even. 
In fact, efc(T) is a sum of determinants of principal submatrices of T in which |A| — k columns and 
corresponding rows are removed. Each such submatrix is naturally associated with a subgraph of A 
— which is necessarily bipartite as well. Therefore, we can conclude that the sign of the determinant 
of every principal submatrix of even order, k, is (— l) fc / 2 . For odd k, such determinants are always 
zero. Warning: The canonical flux distribution need not maximize |efc(T)| for k ^ |A|. The reason 
is that the canonical flux distribution for a subgraph might differ from the one for the full graph 
A. See Theorem 5.1, however. 

(hi). The canonical flux distribution maximizes detT in the bipartite case. It fails, generally, 
to do so in the non-bipartite case; nevertheless, it does have a "maximum property" which is given 
in Theorem A. 2 in the appendix. 

(iv). In the Appendix it is shown that Theorem 3.1 is one of the two key ingredients in a proof 
of Kasteleyn's theorem. 

Proof: By definition, the determinant is a sum over permutations of monomials in the matrix 
elements of T, each of which is a product of the kind e(7r)ti j7r (i) ■ ■ ■ *|A|,tt(|A|) • Here, £(ir) = ±1 

is the signature of it. Using the cycle decomposition of the permutation it we see that the above 

k 

monomial can be written as f\ (— 1)^ _1 11 ^> where Ci, . . . , C k is a family of circuits with the 
property that every vertex of the graph is in precisely one circuit. Here £j denotes the length of 

the circuit Cj. By the definition of the canonical flux distribution, Y[T = Y\ \^xy\ exp[±Z7r/j/2], 

Cj Cj 

where fj is the number of interior triangles of Cj, and the sign in the exponent indicates the 

orientation of Cj. Thus, the determinant is now a sum over all circuit decompositions of terms 

fc 

of the form ] [ J~[ \t xy \ (— 1)^' _1 cos(7r/j/2). Note that the factor 2 is counted by distinguishing 

J'=l Cj 

circuits of different orientations. 
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By Lemma 2.3, /_,- = £ j + 2v j -2. Thus, cos(7r/ 3 -/2) = (-l)^'" 1 cos^lj /2+ttvj -vr) = 

(— l) £j ' cos(ir£j/2 + 7rUj). If is odd, the cosine vanishes. Hence, only even circuits contribute to 
the determinant. This is a crucial property of the canonical flux distribution! Moreover, since every 
vertex must belong to a circuit, and every circuit has even length, vj is also even for all j and hence 
2vj = (mod 4) and does not contribute to the sign of the monomial. Therefore the monomial 

equals ft cos(vr^/2) U M = ft (-1)^ /2 II M = (-1) |A|/2 U M since ELi4 = l A l- N °te 

j = l Cj j=l Cj Cj 

that when |A| is odd there is at least one circuit of odd length in every circuit decomposition, and 
hence detT = 0. 

The last step is to derive relation (3.2), which is geometrically "obvious". It suffices to note in 
our case that D(T) ■ D(T) is a sum of terms, each of which is of the form T>\T>2, where T>\ (likewise 
T> 2 ) denotes a single term in (3.1) corresponding to a single dimer covering. If we superimpose the 
two coverings we get a collection of disjoint circuits C\, . . . , Ck on A. Each site of A is in exactly 
one of these circuits. Additionally, each circuit will have an even length. This "circuit covering" 
of A corresponds to a term in detT. Conversely, each term in detT corresponds to a "circuit 
covering". (Note: it is here that we use the fact that only circuits of even length contribute to 
det T, for otherwise some terms in det T might give rise to "circuit coverings" that contain circuits 
of odd length.) All that is needed is to check that the weights in D(T) 2 correspond to those in 
detT. The weight of a "circuit covering" in detT is 2 n , where n < k is the number of circuits 
whose length exceeds 2. The factor of 2 comes from the two possible orientations of the circuit or, 
in other words, the contribution of a cyclic permutation and its inverse. The same factor 2 n arises 
in D(T) 2 because each circuit can be decomposed into a dimer covering of the circuit in exactly 
two ways. ■ 

IV. RINGS WITH ARBITRARY WEIGHTS 

We begin our study of the problem of maximizing eigenvalue sums of T with respect to fluxes 
by considering the simplest possible case. In the process some notation and identities will be 
established that will prove useful in later sections of this paper. 

A ring of R > 2 vertices (or R edges) is a graph A with |A| = R vertices labeled 1 up to 
R and with edges (1, 2), (2, 3), . . . , (R — 1, R), (R, 1). The hopping matrix is then determined by 
R complex numbers ti2, ■ ■ ■ ,tm with magnitudes given a priori as l^j+il. Note that A is not 
necessarily bipartite, i.e., R = |A| does not have to be even. 

Although the spectrum of T is easy to compute explicitly if is independent of i, and 

hence one might think that our main theorem here, 4.1, is without content, we draw the reader's 
attention to the fact that we shall consider all possible T's. In other words, we shall be dealing with 
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the "random one-dimensional Laplacian" whose spectrum is the object of much current research. 
From this point of view, it is somewhat surprising that some physical quantities of this random 
system can easily be maximized with respect to the flux. 

While our goal is to compute E (T) in (1.2), we shall consider more general functions of the 
eigenvalues of T. Let / : IR + — > IR be a real valued function defined for nonnegative reals, and 
define F by 

|A| 

F(T) = Tr/(r 2 ) = J]/(A J 2 ), (4.1) 

j=i 

where Ai > A2 > ■ . . > A|a| are the eigenvalues of T. The / needed for E is 

fi(x) = (4.2) 

while for In | det T\ it is 

f 2 (x) = \\nx. (4.3) 
Still another physically important function is 

/ 3 (x) = In cosh y^- (4.4) 

appropriate to the free energy, T = —f3~ 1 lnS, in the grand canonical ensemble^" (1.3): 

T = -2/?- 1 Trln(e- /3T + 1) = -2/?- 1 Tr{ln[cosh(/3T/2)] + In 2 - (3T/2) 

= -2/?- 1 Tr/ 3 (/3 2 T 2 /4) -2/3" 1 |A|ln2 (4.5) 

since TrT = 0. Here, (3~ x = (Boltzmann's constant) x (temperature). 

All these functions have the property of being concave, i.e., /(Ax + (1 — \)y) > \f(x) + (1 — 
\)f(y) for all x > 0,y > and < A < 1. (In fact they are strictly concave, i.e., equality implies 
that x = y.) 

These three functions also belong to a more restricted class of functions which we call inte- 
grated Pick functions. These are functions with the integral representation 

00 

f(x) = clnx + J In (l + fi(ds) (4.6) 



t In (4.5) we have set the chemical potential /i equal to zero. For a bipartite lattice this yields an average particle 
number M = 2N = |A|, which follows from M = 2 exp(-/3\j) [1 + exp(-(3Xj )] _1 together with the 

(A, — A) pairing. 
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where c > and where \i is a nonnegative measure on [0, oo) such that this integral is finite. For 
the functions in (4.2)-(4.4) we have this integral representation with \i{ds) = (const. )s -1 / 2 ds and 
c = for (4.2), n(ds) = and c = 1 for (4.3) and /i(ds) = 

T,kLo 6 ( s ~ + l)]" 2 ) ds > and c = with 5 = Dirac's 5-function. See [KL, eqs. (3.12), (3.16)]. 
Combining (4.1) with (4.6) yields 

oo 

F(T) = clndetT 2 + J In det(l + T 2 /s)n(ds), (4.7) 
o 

and we see that the problem of maximizing F(T) is reduced to that of maximizing various deter- 
minants with respect to the flux. 

The function G(T) given in (1.4) cannot be represented in the form (4.7); nevertheless we 
shall also be able to maximize G(T). 

4.1. THEOREM (Maximizing flux for the ring). Consider a hopping matrix T with 
arbitrary, but fixed amplitudes \t xy \ on a ring of R sites, let f be an integrated Pick function given 
by (4-6) and let F(T) be as in (4-1)- Then the canonical flux 

tt(R + 2)/2 (mod27r) maximizes both F(T) and the gap G{T) if R is even. If R is odd, F{T) and 
G{T) are maximized by both of the choices and tt. 

Remark: When R is odd the canonical flux is always tt/2 or 37r/2 and never or tt. 

Proof: For F(T) it suffices, by formula (4.7), to show that the flux described above maximizes 
det(c 2 + T 2 ) = | det(zc + T)| 2 for all real numbers c. 

As a first step we observe that for the invariants (or elementary symmetric functions) we have 
that ejfe(T) = for k odd and 1 < k < R— 1. This follows directly from remark (ii) after Theorem 
3.1 when R is even. When R is odd it also follows from remark (ii) together with the observation 
that every proper subgraph of a ring is bipartite. We also see, from remark (ii), that the sign of 
e 2m (T) is (-l) m . 

With this information about the signs of the e^'s, we can write, from (3.4) with z = ic, 



det(c 2 



l A l/ 2 , A , \ 

E |e 2m (T)|cl A l- 2m , |A| even 

m=0 / , , . 

) s 2 (4.8) 

'(|A|-l)/2 V 

E |e 2m (T)|cl A l- 2m +(detT) 2 , |A|odd. 

m=0 / 



For future use we remark that both parts of (4.8), hold for any bipartite graph A, not just a ring 
with an even number of sites. 
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As a second step we shall show that in the expression (3.4) for det(T + z), with z ed, the 
invariants e^(T) are independent of the flux if k < R. This is true only for a ring. Note that the 
invariants are real since T is Hermitian. Recall that the number efc(T) can be computed from T 
by calculating the subdeterminants of T with any R — k columns and corresponding rows removed, 
and then summing these numbers over all possible removals. The result is a sum of monomials of 
the form n( — 1)^ _1 T\T where the product is taken over all circuits Cj that cover the subgraph 

obtained by removing k vertices and the corresponding edges. But for k < R the only circuits that 
cover this subgraph form a dimer covering; their contribution does not depend on the flux but only 
on the numbers l^j+il. Thus the only term in (3.4) that depends on the flux is e^\(T) = detT. 

In both cases in (4.8), the problem of maximizing det(c 2 + T 2 ) is seen to be the same as 
maximizing | detT|. If R is even this problem is solved in Theorem 3.1. If R is odd there are 
precisely two circuits that contribute to det T. These are the circuits that traverse the entire ring 
(in either direction) and correspond to an even permutation. Thus, for a ring of odd length 



from which we see that <3? = or $ = ir maximizes | detT|, and hence also F(T). This completes 
the proof for F(T). 

To compute G(T) we return to (4.7) and write Q$(A) = det(T - A) = P{\) + detT, with 
P being a polynomial of order R whose coefficients are independent of the flux P is even if 
R is even and P is odd if R is odd. We note that as detT varies between its maximum and 
minimum values, Q$ always has R roots. We leave it to the reader to verify the following with 
the aid of a graph of Q$(A). Even R: The maximum separation between Xr/2 and A#/ 2 +i is 
achieved by making |Q$(0)| as large as possible. Since Q$(0) = detT, this means choosing <3? 
to make |detT| as large as possible — as stated in our theorem. Odd R: If detT = 0, the 
eigenvalues of Q$ are paired (because P is odd). Then —X(r-i)/2 = X(R +3 y 2 and \(r+i)/ 2 = 0. 
Thus G + := [\( R+ i)/2 - \r+s)/2] = [\r-i)/2 ~ \r+i)/2] =■ G- when detT = 0. As |detT| 
increases, either G + increases, G_ decreases and \(r+i)/2 > or vice versa and \(r+\)/2 < 0. 
Thus G = max(G_)_,G_) and this increases with | detT|. By (4.9) we see that <1> = or <I> = ir 
maximizes | det T\. ■ 

V. TREES OF RINGS 

Most of the results in Theorem 4.1 for bipartite rings with arbitrary hopping amplitudes \t xy \ 
can be extended to a much larger class of planar graphs. Two special cases of this class are the 




(4.9) 



16 



EHLML-21/Sept /92 



ladders and the necklaces; they are discussed in detail in the next section because even stronger 
results can be obtained for them. It was those two classes, in fact, that were the origin of this work 
and that were reported in [LE]. 

A planar graph A is said to be a tree of rings if and only if A has an embedding in 1R 2 such 
that every circuit in A has no interior vertices. 

The simplest example consists of two rings which have exactly one vertex in common. Another 
example consists of two rings that have exactly one edge (i.e., two neighboring vertices) in common. 
More generally, one can have a "tree of rings" in which two successive rings share either one edge or 
one vertex. The canonical flux distribution for a tree of rings would have flux (tt/2)[(£ — 2)(mod4)] 
in each circuit of length £. 

5.1. THEOREM (Maximizing flux for bipartite trees of rings). Let A be a bipartite, 
planar graph that is a tree of rings and let \t xy \ be arbitrary given hopping amplitudes. For f an 
integrated Pick function, let F(T) be as in (4-1)- Then the canonical flux distribution maximizes 
F(T). Moreover, it also maximizes the magnitude of each elementary symmetric function efc(T) 
defined in (3.3). 

Proof: As in the proof of Theorem 4.1 (c/. eq. (4.8)) we have that det(c 2 + T 2 ) will be 
maximized if we can simultaneously maximize all the |e/c(T)|'s and if they all have the sign (— l) fc / 2 . 
The latter question was dealt with in the remark (ii) just after Theorem 3.1. 

Each efe(T) can be evaluated as a sum of determinants of principal submatrices of T of order 
k. In terms of graphs, a particular term in the sum is the determinant of T restricted to a subgraph 
A' with |A'| = k. The important point is that the circuits of A' are (i) a subset of the circuits of 
A and (ii') they have no interior points. The canonical flux distribution for A' is the same as for 
A; this means that if C is a circuit that is both in A and in A' then &c = &c = or 7r where $c 
is the canonical flux through C (in A) and <&' c is the canonical flux (in A'). (Note: The only way 
in which $<7 could differ from <&' c is if C had some interior vertices that were removed in passing 
from A to A'. But C had no interior vertices to start with.) Hence each subdeterminant appearing 
in efc(T) is maximized (in absolute value) by the original canonical flux distribution in A. Since 
the signs of all these subdeterminants are the same, in fact they depend only on k (see remark (ii) 
after Theorem 3.1), we see that |efc(T)| is maximized. ■ 

VI. LADDERS AND NECKLACES 

Most, but not all the graphs considered in this section are special cases of those discussed in 
Section 5. Here we consider certain graphs that are finite subsets of the infinite lattice Z 2 , which 
is the infinite embedded graph whose vertices are points in the plane with integer coordinates and 
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whose edges are the horizontal and vertical line segments joining vertices a unit distance apart. 
Of particular importance are boxes, which are the subgraphs of Z 2 with 4 vertices and 4 edges 
forming a circuit. In general, our graphs need not be subgraphs of Z 2 , i.e., they need not contain 
all the edges of Z 2 that connect the vertices V in our graph. (Example: A contains the 4 vertices 
of a box but only 3 of its edges.) Evidently, all our graphs are bipartite — with the A — B 
decomposition of their vertices being the one inherited from Z 2 . Before describing them in detail, 
a few remarks are needed. 

If T is a hopping matrix of a bipartite graph A, the matrix T 2 is evidently block diagonal, 
i.e. T 2 = ( a ^ „ ), where ax is \A\ x \A\, fix is \B\ x \B\ and both are positive semidefinite. 



fi T 

Assuming that v := \B\ — \A\ > 0, we have that the eigenvalues satisfy 

spec(fix) = spec(aT) U {y zeros}. (6.1) 
This is a simple consequence of the fact that T = ( ® ) , so that ax = MM* and (3t = 



^M* / 

M*M. Since the eigenvalues of T come in pairs, and there are |A| of them, we conclude that 

spec(T 2 ) = {spec(aT) with double multiplicity} U {v zeros}, (6.2) 
spec(T) = spec (\/ax) U spec {—\fax) U {y zeros}. (6.3) 

Thus spec(T) is determined by either ax or fix alone. 

The matrix ax has diagonal elements, (ax) aa = J2btB \tab\ 2 for a & A. The off-diagonal 
elements of ax can be thought of as a hopping matrix of a new graph Aa, which need not be 
planar. The vertices of A^ are the yl-vertices of A. A necessary condition for the pair (a, a') 
to be an edge of A^ is that there is a .B-vertex, b, such that (a, b) and (b, a') are edges in A. 
There may be more than one such b for a given pair (a, a') and it is important to note that since 
(ax)aa' = ^2i,tabtba' , it c & n happen that (ax)aa' = 0, in which case (a, a') is not an edge of A^, 
in conformity with our earlier convention. Similar remarks hold for {(3x)bb' and A#. 

There are special edges in A^ or A# which we call interior diagonals. These are edges (a, a') 
or (6, b') in which a and a 1 (or b and b') belong to some (same) box S that is a subgraph of A. This 
set of edges is denoted by Da- Since circuits of length 4 can only be the edges of boxes, it follows 
that the interior diagonals are the only edges that can possibly disappear from A^ because of the 
equality ^2 b t a btb a > = 0. Likewise for A^ and Db- 

The graphs we shall consider here can be described as follows. Let A^ = A^ ~ Da, i.e., 
the vertices of A' A are those of A^ but the edges are those of A^ without the interior diagonals. 
Analogously, A' B is defined. We say that A is a hidden tree if either A^ or A' B is a tree (i.e., does 
not contain any circuits). 
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Two important examples introduced in [LE] are ladders and necklaces. Each is a connected 
union of n boxes, labelled 1,2, ...,n, forming a one-dimensional array. In the ladder the boxes 
are joined along (parallel) edges, with box j connected to j + 1. In the necklace, boxes j and 
j + 1 have only a single vertex in common, and these vertices are either all A or all B. In both 
examples, boxes j and k are disjoint if \j — k\ > 1. With the usual orientation of Z 2 , ladders are 
either horizontal or vertical, while a necklace runs at 45° to either of these directions. 

One can generalize the ladder by allowing occasional 90° bends, while still keeping the one- 
dimensional character. Now, squares j and k can now have a vertex (but not an edge) in common 
if \j — k\ = 2; squares 1 and n are disjoint. The bends cannot be completely arbitrary, however, 
because the "hidden tree" condition must be maintained. It can never be maintained for a necklace 
with 90° bends. 

Another, not so trivial example is that in which A is the union of 4 squares, all of which 
have a vertex in common and which together form a square of side length 2. Here \A\ = 4 and 
\B\ = 5. Although A^ is a tree, A^ is not a tree. (6.1) notwithstanding, it is somewhat surprising, 
when viewing the graphs for A a and A# , that they have the same spectrum — except for one zero 
eigenvalue. Another example in this vein is the A that resembles a (3,2) Young diagram, i.e. 3 
squares in a horizontal row and 2 squares, also in a horizontal row, directly beneath them. The 
simplest case that is not a hidden tree, and for which none of our theorems apply, is two rows of 
three squares each. 

Finally, we make some remarks about the next theorem. 

(i) . Ladders and necklaces are trees of rings, but the last three examples (two rows of squares) 
are not. Thus, for ladders and necklaces, the fact that the canonical flux distribution maximizes 
Tr|T|, T and |detT| is already covered by Theorem 5.1. The statement about the gap is new, 
however, as is the method of proof. 

(ii) . The concave function F{T 2 ) is definitely a generalization of the function in (4.1). Not 
all concave functions (even those that are invariant under unitary transformations) are eigenvalue 
sums as in (4.1). In particular, the sum of the k lowest eigenvalues of T 2 is not such a function 
and it is needed, in fact, to prove the theorem about the gap. 

6.1. THEOREM (Canonical flux maximizes concave functions on hidden trees). 
Let A be a graph that is a subset of Z 2 and suppose that A' A (resp. A' B ) is a tree. Let F be a 
concave function on the cone of positive definite matrices of order \A\ (resp. \B\) with the property 
that F(U*PU) = F(P) for every P > and every gauge transformation U (restricted to Aa, of 
course). Finally, let {\t xy \} be unit hopping amplitudes on A, i.e. \t xy \ = 1 if (x,y) 6 E(A). 
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Our conclusion is that among all hopping matrices T with this unit hopping amplitude, F(cxt) 
(resp. F(Pt)) is maximized by the canonical flux. 

Proof: We shall assume A' A is a tree. The proof for A' B is similar. Assume T maximizes F(ckt)- 
Let U = {u x 6 xy } XtyeV (A) be the following gauge transformation: If x = (n,m) with n, m 6 Z, then 
u x = (—1)™. Let Y = U*TU. By concavity and gauge invariance, the block diagonal matrix 
P 2 := \{T 2 + Y 2 ) satisfies F(a P ) > \F(a T ) + \F(a Y ) = F(a T ). This inequality proves our 
theorem if we can show that the matrix ap can be achieved by the canonical flux distribution, i.e. 
if there is a gauge such that C := (T with the canonical flux distribution) satisfies «c = ap. 

Now note that T 2 and Y 2 = U*T 2 U are related as follows: 

(Y 2 )xy = —(T 2 )xy if (x,y) is an interior diagonal 
{Y 2 ) xy = (T 2 ) xy otherwise. 

Therefore, (P 2 ) xy = for (x,y) G D A U D B and {P 2 ) xy = (T 2 ) xy otherwise. 

For any gauge, (C 2 ) xy = if (x, y) is an interior diagonal. This is so because the flux through 
each box in A is ir, and if we label the four vertices 1, 2, 3, 4 (counterclockwise) with 1 and 3 6 Va 
we have (C 2 )i 3 = Ci 2 C 2 -i+CuC^. But C12C23C34C41 = -1 and C14 = C41 = I/C41, C43 = 1/C 34 
(since, e.g., | C14 1 2 = 1), so (C 2 )i3 = 0, as required. As for the diagonal elements, they are clearly 
equal, i.e., (C 2 ) xx = (P 2 ) xx . 

Finally we have to compare the other matrix elements of C 2 and P 2 on A^. In fact, they 
are both nonzero only on A^, in which case they satisfy |(C 2 ) aa /| = |(-P 2 ) a a'| = 1 because there is 
precisely one path (i.e., .B-vertex) between a and a'. Therefore (C 2 ) aa > = exp[i0(a, a')](P 2 ) aa > for 
each edge in A^' • The relevant question is then the following: Is there a gauge transformation U 
such that {U*C 2 U) aa ' = (P 2 )aa'^ In other words, can we find u x = exp[i(f) x ] such that <j) a — (f) a > = 
6(a,a') for every edge (a, a') in A' A ? Since A' A is a tree, the answer is trivially, yes. All one-forms 
on a tree are exact. ■ 

Applications: For a graph with hopping amplitudes satisfying the hypotheses of Theorem 6.1 
the canonical flux distribution yields: 

(a) The lowest ground state energy E (T) and free energy F(T) for all temperatures. 

(b) The largest | detT| and gap G(T). 

The functions —Eo(T), —T(T) and log | detT| are concave since they are integrated Pick functions, 

which are concave in T 2 as mentioned in Section IV. The gap G(T) can be computed from the 

matrices ay and (3t- Assuming that v = \B\ — \A\ > we see that G(T) = 2(infspec ar) 1 ' 2 and 
fu+l \V 2 

that G{T) = 2 I Yl 7j J , where 7, denote the eigenvalues of (5t arranged in increasing order. Of 
course we used that 71 = 72 = • • • = 7^ = 0. Now the sum of the first k eigenvalues of a Hermitian 
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matrix H is a concave function of H. Moreover since x i-> y/x is concave and increasing, we see 
that G(T) is a concave function of «t (resp. (3t)- Moreover it is gauge invariant and hence satisfies 
the assumptions of Theorem 6.1. Note that we had to discuss the above two formulas for G(T) 
since both possibilities (A^ is a tree or A' B is a tree) have to be considered. 

In the case of ladders and necklaces the result about _E (r) , JF(T) and | detT| were covered 
by Theorem 5.1, but the examples with two rows of boxes, cited above, were not covered. (In 
fact, the two-rowed examples cannot be extended to the full generality of Theorem 5.1; see Section 
VIIB.) For ladders and necklaces, the result about the gap is not covered by Theorem 5.1. In [LE] 
it was mistakenly asserted in Theorem 1 that (a) and (b) hold for fully generalized ladders and 
necklaces in which arbitrary 90° bends are allowed. Indeed, for Eq(T),!F{T) and |detT| this is 
correct (by Theorem 5.1). For the gap, however, we must use Theorem 6.1 and this fails for bent 
necklaces and for ladders with arbitrary bends. It does hold, however, for generalized ladders that 
are hidden trees. 

Additional remarks and examples: There are two more cases where the concavity argument 
in the proof of Theorem 6.1 is applicable but where the graph is not necessarily a hidden tree, or 
even planar. 

(i) Row of cubes: Instead of a row of squares as in the ladder, take A to be a row of cubes 
joined on their faces, i.e., neighboring cubes have four edges and one face in common. Such a graph 
is, in fact, bipartite and planar, but it is not a hidden tree. If there are n cubes then A is the 4 x n 
planar, square lattice with "periodic boundary conditions" in one direction. [Indeed, we can even 
make the row of cubes into a torus (i.e., attach the first cube to the last), which is the same thing 
as the 4 x (n + 1) planar, square lattice with periodic boundary conditions in both directions; the 
following argument will continue to work in this case provided n is odd.] 

We assume, as in Theorem 6.1, that \t xy \ = 1 for every edge. The flux in every face can easily 
be arranged to be tt in the following way. Start with the face that cube 1 and cube 2 have in 
common and put flux tt through it by making t xy = 1 on three edges and —1 on the fourth. Then 
use the negative of this on the corresponding edges that cube 2 and cube 3 have in common - 
and so on alternately. Finally, set t xy = +1 on the remaining edges, i.e., those edges that are 
perpendicular to the faces between the cubes. A first application of the concavity argument shows 
that we get an upper bound in terms of T 2 , but without interior diagonals on the faces common to 
all the cubes. In a second, similar application of the argument we can get an upper bound in terms 
of a matrix that has no interior diagonals on any of the other faces of the cubes as well. In fact the 
only nonzero elements of T 2 that remain will consist of four independent one-dimensional chains 
(or else, in the case of the torus, four independent rings of length n + 1 with zero flux through 
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each ring). Theorem 4.1 says that if n + 1 = 2 (mod 4) the optimizing flux for such rings is zero 
and hence the above concavity argument shows that our choice of fluxes cannot be improved. If 
n + 1 = 0(mod 4) the optimum choice for such a ring is flux tt. This can be achieved by a slight 
modification of our initial choice of the t xy s along the edges perpendicular to the inter-cube faces. 
Initially we chose them all to be +1 but now we choose them all to be +1 except for the n th cube. 
There we choose t xy = — 1 along the four perpendicular edges. 

(ii) SU (2) -valued fields: We have considered the case that t xy = \t xy \ ex.p[i6(x,y)] with e l8 the 
unknown variable. It is also amusing to replace e l6 , which is in U{1) by a 2 x 2 matrix U xy in SU(2). 
In other words, T becomes a 2|A| x 2|A| Hermitian matrix in which each t xy (for x,y £ ^(A)) 
equals a given number \t xy \ times an (x, y)-dependent element of SU(2). We require t xy = t* y . 
Theorem 6.1 goes through in this case when \t xy \ = 1. We do not know whether Theorem 3.1, for 
instance, can be generalized to the SU{2) case. 

There is, however, an interesting special feature of SU(2). For the ladder or row of cubes, we 
had to break the translation invariance from period one to period two in order to achieve flux tt in 
each face, i.e., T could not be made the same in every box but we had to translate by two boxes 
(or cubes) in order to recover T. With SU(2) fields we can achieve the optimal flux distribution 
with period one. This means the following. We require that the product of the four t xy s around 
a square face (which is now a matrix product, of course) is the matrix —I G SU{2). This can be 
achieved by placing ia 1 along all horizontal edges, ia 2 along all vertical edges and (in the case of 
cubes) ia 3 along all the edges in the remaining direction. Here a 1 , a 2 , a 3 are the Pauli matrices. 
We then have —I in every face because, for example (ia 1 )(ia 2 )(—ia 1 )(—ia 2 ) = —I. 

VII. SOME CONJECTURES AND COUNTEREXAMPLES. 

(A) . The smallest determinant: A natural question, to be compared with Theorem 3.1, 
is "Which flux distribution minimizes |detT| for a bipartite lattice?" Since the canonical flux 
distribution maximizes | det T\ and since it places flux tt in each square face (which is the maximum 
possible flux), it might be supposed that the answer to the question is zero flux, i.e. set t xy = \t xy \ 
for every edge in E. In the case where A is a simply connected net of boxes on Z 2 and \t xy \ = 1 
the determinant was computed by Deift and Tomei [DT] to have the three possible values 0,-1 or 
+1. Despite this supportive example the above conjecture is wrong. In the case of two boxes with 
one common edge and \t xy \ = 1, the determinant vanishes when the flux in each square is tt/3. On 
the other hand, detT = — 1 when the flux is zero. 

(B) The smallest energy: In Section V we have seen examples of some graphs whose energy 
is minimized by the canonical flux distribution for arbitrary hopping amplitudes. Moreover, | det T| 
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is always maximized by that flux distribution for every bipartite, planar graph with \A\ = \B\. For 
such an arbitrary graph it is therefore natural to conjecture that E(T) is also always minimized 
by the canonical flux distribution. Alas, this conjecture is also false for arbitrary \t xy \, as we how 
show by an example. 

In Z 2 consider the graph consisting of four boxes arranged in a square, i.e. A has the nine 
vertices o = (0, 0), a = (1, 0), b = (0, 1), c = (-1, 0), d = (0, -1) and (±1, ±1). If the conjecture 
were true for arbitrary amplitudes then E (T) would always be minimal if the flux in each square 
is 7r. If we now let \t oa \, \t b\, \t c\ an d \t d\ tend to zero, E (T) becomes that of a ring of eight 
sites for which Eq(T) is minimized by flux tt (not = 4tt), as we saw in Theorem 4.1. Indeed, to 
see that and tt do not give the same value for Eq(T) in general, for this ring, assume that the 
i's on the ring all have amplitudes equal to one. If the flux is zero then Eq(T) = — Tr\T\ is easily 
computed to be —2(1 + \f~2) and for flux tt it equals — 2 5 / 4 (-\A + V% + V V% — 1) which is more 
negative. 

VIII. THE FALICOV-KIMBALL MODEL 

So far all our models concerned hopping matrices on a graph A. It is likely that the results 
obtained do not hold if we add a diagonal term to T, i.e., if we replace T by the matrix 

H = -T + 2UW. (8.1) 

Here W is a real diagonal matrix satisfying < W x < 1 for each x and U is a given real number 
called the coupling constant. 

The eigenvalues of H can be interpreted as the possible energy levels of a single electron 
hopping on a graph A with kinetic energy — T and potential energy 2UW. This model, with W x 
restricted to be or 1, was introduced in [FK] as a model for a semiconductor- metal transition, 
and it was studied extensively in [KL], where it was called the "static model", and in [BS]. Our 
generalization to W x 6 [0, 1] for each x is a mild one that we include here primarily because it 
can be handled without extra complication. The points of view taken in [KL] were different from 
that in [FK]. There, the model was considered either as a simplified version of the Hubbard model 
in which electrons of one sign of spin are infinitely massive and therefore do not hop, or else as a 
model of independent electrons interacting with static nuclei. In the first view, U > and U < 
are both relevant. In the second, U < is the physically relevant sign, and W x = 1 (resp. 0) 
denotes the presence (resp. absence) of a nucleus at x. 
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The eigenvalues of H in (8.1) are denoted by —v (not —A) which, as usual, are ordered 
v\ > v-i > . . . > z/|A|. The ground state energy of a system of ./V spinless electrons, interacting with 
a magnetic field and with the nuclei, is given, as usual, by 

N 

4"> = -J> (8.2) 

i=l 

(There is no 2 here, as in (1.1), because there is no spin.) Notice that the spectrum of (8.1) is still 
invariant under gauge transformations. Thus the energy Eq N ^ depends on N, the flux distribution 
and, of course on W. 

As mentioned before, minimizing the energy over all fluxes with a fixed W is presumably a 
hopeless endeavor but the situation becomes easier if one tries to minimize the energy with respect 
to the fluxes and W. It is clear that the minimum is attained when W = if U > or when W = I 
if U < 0, which is uninteresting both mathematically and physically. If, however, we introduce 
N n := Yl W x = (total charge of the static particles), then the half-filled band condition (from the 
Hubbard model point of view — at least) is that N n + N = |A|. This is the case that parallels the 
restriction in the previous parts of this paper, since it means that on the average each lattice site 
is occupied by one particle. In fact we shall be a bit more general in the U > case, and will treat 
a slightly different case when U < 0. We shall optimize the energy over all flux distributions and 
all potentials W and all choices of ./V subject to one of the following three constraints. 

J2W X + N<2\A\ HU<0 (8.3a) 
J2 w x + N <2\B\ HU<0 (8.36) 

x£A 

J2 W * + N > l A l if C/ > 0. (8.3c) 

x£A 

Theorem 2.1 in [KL] says that in these three cases we can easily compute the minimum of 
Eq N ^ with respect to W and M — regardless of the flux distribution. This result requires only 
one particular structure of the graph, namely that it is bipartite. Nothing else is required. As we 
shall see in Theorem 8.2 the result is that the nuclei want to occupy only the A sites or the B 
sites in order to minimize the total energy. We emphasie again that this fact is independent of the 
magnetic field. 

Lemma 8.1 and Theorem 8.2 are really a transcription to the W x € [0, 1] case of Lemma 2.2 
and Theorem 2.1 in [KL]. Since they are short we give them here. It is convenient to introduce the 
matrix S = 2W - I, so S is diagonal with S x G [-1,1]. Thus H = h + UI with 

h = -T+US. (8.4) 
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The matrix h has eigenvalues —fii < —fi2 < ■ ■ ■ < — A*|A| 5 with = za,- + £/. 

8.1. LEMMA (Maximization with respect to the nuclear configuration). Let A be a 

bipartite (not necessarily planar) graph and T a prescribed hopping matrix on A. We consider all 
functions S on the vertices of A satisfying — 1 < S x < 1 for all x. Let F be a concave, nondecreasing 
function from the set of Hermitian positive semidefinite matrices into the reals. Assume also that 
F is gauge invariant, F(U*PU) = F(P) when U is a gauge transformation. 

Then F{h 2 ) is maximized with respect to S at S = V and at S = —V, where 

(+lforxeA 
Vx ~\-l forxGB ■ { *- b) 

Lf F is strictly concave and strictly increasing then these are the only maximizers. 

Proof: The matrix V = V* is a gauge transformation and hence h and h! := VhV satisfy 
F(h 2 ) = F(h' 2 ). Now h 2 = T 2 + U 2 S 2 - U(TS + ST) and, since VTV = -T and VSV = S, 
h' 2 = T 2 + U 2 S 2 + U(TS + ST). By concavity 

F(h 2 ) = \ [F{h 2 ) + F(h' 2 } < F(T 2 + U 2 S 2 ) < F(T 2 + U 2 L), (8.6) 

since F is nondecreasing and S 2 < L. Note that T 2 + U 2 I = h 2 when S is chosen to be +V or —V. 
If F is strictly increasing and strictly concave we can have equality in (8.6) only if TS + ST = 
and S 2 = L. The former implies that t xy (S x + S y ) = 0, which implies (since A is connected) that 
S = (constant)y. The latter implies that (constant) = ±1. I 

8.2. THEOREM (Energy minima with respect to nuclear configurations). Let A be 
a bipartite graph (not necessarily planar) and let T be a prescribed hopping matrix. We consider 
functions W on the vertices of A satisfying < W x < 1. 

For the three cases given in (8.3) the minimum value of with respect to W and N is 

uniquely achieved as follows 

N = \A\ and W = W A := \{I + V), U < (8.7a) 

N = \B\ and W = W B := \{I - V), U<0 (8.7b) 

N = \B\ and W = W A or N = \A\ and W = W B , U > 0. (8.7c) 

For each of these three cases, the minimum E^ is given by 

E { N) = -±Tr\h\ + ±Tih + UN, (8.8) 

where h = -T + U{2W - L). 
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Proof: By (8.2) and (8.4) we have that E = - J2f=i v 3 > ~^\h\ + ^Tvh + UN =: A(h). By 
Lemma 8.1, Tr|/i| = TrVtv 2 is maximal precisely at W = Wa or W = Wb since x — > is a strictly 
increasing and strictly concave function. Also, ±Tr/i + [/A = U{£ x W x - |A|/2 + A} =: If 
17 >0, B(h) > |A|/2. If C7 < B(h) > U(2\A\ - |A|/2) for (8.3a) and B(h) > U(2\B\ - |A|/2) for 
(8.3b). These three lower bounds in B(h) are attained (under conditions (8.3)) if A and W satisfy 
(8.7). 

To complete the proof, we have to show that the lower bound on — ^Tr|/i| =: C(h), given 
in Lemma 8.1, is compatible with the condition on A given in (8.7) (when W is also that given 
in (8.7)). For example, we have to show that if W = Wa and U < then the sum of the 
negative eigenvalues of h (namely — >o A*j = — ^ Tr | /i | + ^Tr/i) equals the sum of the lowest \A\ 
eigenvalues of h. In other words, we have to show that h has exactly \A\ negative eigenvalues. We 
do so now with a proof different from the one in [KL]. First note that for t G [0, l],h t = —tT + UV 
has no zero eigenvalues when U ^ since 

h 2 = t 2 T 2 + JJ2J > JJ2J > _ Second, the matrix 
h = UV has precisely \A\ negative eigenvalues because U < 0. Since the eigenvalues of h t are 
continuous functions of t and because no eigenvalue can cross zero, h also has precisely \ A\ negative 
eigenvalues. The other two cases (8.7b) and (8.7c) are treated in the same fashion. ■ 

The next theorem is our main result about the FK model in a magnetic field. 

8.3. THEOREM (Canonical flux minimizes energy on trees of rings). Let A be 
a bipartite tree of rings, T a hopping matrix with arbitrarily prescribed amplitudes, \t xy \. As in 
Theorem 8. 2 we consider functions W on the vertices of A satisfying < W x < 1 . 

For the three cases given in (8.3) the minimum value of with respect to the flux distri- 

bution, N and W is achieved by the canonical flux distribution together with the N and W given 
by (8.7). 

Proof: Minimizing first with respect to N and W, we can assume, by Theorem 8.2, that (8.7) 
is satisfied. In these cases we have that E = — ^Tr\fiA,B\+ constant, where the constant depends 
on the case but not on the flux distribution. In each case, W = Wa or Wb and we denote the 
two choices of h by Ka and hs- Note that N no longer enters the discussion. Our only goal 
now is to maximize Tr|/i^ )B | with respect to the flux distribution. But Tr|/i^ !B | = Tr^Jh A B and 
h AB = T 2 + U 2 I since S AB = I. The function x i— > yfx is an integrated Pick function, i.e., 

oo 

\fx = d j ln(l + xj ' s)s~ 1 / 2 ds for some constant d > (see. 4.6). Hence maximizing Tr|/iA,s| is 
o 

reduced to maximizing det(c 2 + h 2 A B ) = det(c 2 + U 2 +T 2 ) for all constants c. That this is achieved 
by the canonical flux distribution on trees of rings is precisely the content of Theorem 5.1. ■ 
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APPENDIX: KASTELEYN'S THEOREM 

We give here a different, and we believe more transparent proof of a deep theorem due to 
Kasteleyn [KP], which is one of the main tools for counting dimer configurations on planar graphs. 
Let us emphasize that A is now a finite graph that is not necessarily bipartite. 

Historically, the motivation behind Kastelyn's theorem was an attempt to calculate efficiently 
the partition function D(T) in (3.1) for large planar graphs — by reducing the problem to the 
calculation of a determinant. This was accomplished by Temperley and Fisher [TF] in special 
cases, but independently and in full generality by Kasteleyn [KP]. The starting point was Pfaff's 
theorem for an antisymmetric matrix A (of even order): det A = Pf(^l) 2 . Here Pf(-A) is the 
Pfaffian of A (more precisely, the Pfaffian of the upper triangular array of A = {a X y}i<x<y<N), 
defined by 



where the sum is over all permutations ir 6 5V with 7r(l) < 7r(3) < 7r(5) < . . . and with < 
n(i + 1) for odd i. Also, e(tt) is the signature of ir. 

Each term in Pf(^4) corresponds to a dimer covering of the graph A with N vertices and with 
edges corresponding to the nonzero elements of A. 

The Kasteleyn, Temperley-Fisher idea is to set a xy = exp[i9(x,y)]\t xy \ for x < y, with 9(x,y) 
chosen so that all terms in (A.l) have a common argument 9. Then, trivially, D(T) = e 1 Pf(-A) 
for some 9 and | det A\ = D(T) 2 . 

In Theorem 3.1 we solved the D(T) problem, from a different perspective, by using the canon- 
ical flux distribution: |detT| = D(T) 2 . Moreover, we gave a very simple rule (in the proof of 
Theorem 3.1 and in the remark following Lemma 2.2) for an explicit construction of exp[i(9(x, y)\. 
We did not try to construct an antisymmetric matrix, but it is a fact that among the gauge equiv- 
alent T's with canonical flux distribution, there is one that is antisymmetric. This is Theorem A.l 
below. We note here that detT is a gauge invariant quantity, but Pf(T) is not. 

Theorem A.l, together with Theorem 3.1 yields Kasteleyn's theorem, namely that there is 
always a real, antisymmetric matrix A such that T = iA and Pf(A) = D{T). To see this impli- 
cation, note that Theorem (3.1) says det ,4 = det(tT) = (-l)l A l/ 2 (-l)l A l/ 2 L>(T) 2 . On the other 
hand det A = Pf (A) 2 . We can then make Pf(-A) = +D(T) by multiplying the first row and column 
of A by —1, if necessary. One way in which our proof is a little simpler than Kasteleyn's is that 
graphs with cut-points do not require special treatment, in either of Theorems A.l or 3.1. 

Theorem A. 2 is another corollary of Theorem A.l. It answers the question: For which class of 
matrices (in the non-bipartite case) does the canonical flux distribution maximize | detT|? Clearly 




(Al) 
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the canonical flux cannot maximize | det T\ in general. (A simple counterexample is provided by 
the triangle, |A| = 3 where | detT| is maximized by flux or it (Theorem 4.1) but the canonical 
flux is 7r/2.) 

A.l. THEOREM (Canonical flux and antisymmetry). Let A be a finite planar graph 
with given hopping amplitudes \t xy \. There exists a gauge such that the Hermitian hopping matrix 
T defined by the canonical flux distribution has purely imaginary elements, i.e., T = iA with 
A T = —A and A real. 

Proof: We can assume A is triangulated and we can start with the remark following Lemma 
2.2 which states (indeed, gives an explicit rule) that there is a matrix T° whose fluxes are canonical 
and such that t° y 6 {i, —i, 1, —1} for all (x, y) € E. Edges for which t° y = ±i will be called "good" 
and those for which t xy = ±1 will be called "bad". Our goal is to find a gauge transformation 
U xy = u x 8 xy with u x €E {i, 1} for all x 6 V, such that U*TU has no bad edges. We shall do so 
by showing that if T is any matrix with canonical fluxes and with t xy £ {i, —i, 1,-1} and with at 
least one bad edge then we can find a gauge transformation U with u x € {i, 1} such that U*TU 
has at least one less bad edge than T has. Since the number of edges of A is finite, the theorem is 
proved by induction. 

Let (a, b) G E be a bad edge. Consider the set of sites S = {y € V : there is a path from a to 
y whose edges are all good }; by definition a 6 S. Set u x = i if x G S and u x = 1 if x g" S. Clearly, 
if (c, d) G E was a good edge for T it remains a good edge for U*TU (because either u c = ua = 1 
or Uc = u d = i). 

To complete the proof we have only to show that the edge (a, b) has become a good edge. 
Since a 6 S, we have to show that b S. Indeed, suppose b 6 S. Then there is a circuit 
C = a, x\, X2, ■ ■ ■ , x n b, of length greater than 2, such that the edges (a, xi), (x\, X2) • • • , (x n -i,x n ) 
are good while (x n ,b) is bad. We claim that this is impossible; in fact every circuit must have an 
even number of bad edges. To see this, use eq. (2.2) in Lemma 2.3, which says that / = I (mod2). 
Here, / is the number of (triangular) faces inside C. We have <E>c = fhrx through C = ±7r//2 (by 
the definition of the canonical flux distribution). On the other hand, &c = ( 7r /2) J2g(^)- Here, 
the sum is over the good edges and +1 or —1 is taken according to the direction in which the edge 
is traversed when C is traversed in an anticlockwise sense. In any event, &c = ( 7r /2){|G'|(mod2)} 
where |G| is the number of good edges in C. Thus, / = |G|(mod2), which proves our assertion. ■ 

A. 2. THEOREM (Canonical flux maximizes antisymmetric determinants). Let A 
be a planar (not necessarily bipartite) graph with hopping amplitudes \t xy \ given. The canonical 
flux distribution maximizes \ detT| among all flux distributions such that T is both Hermitian and 
antisymmetric. 
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Proof: If T is antisymmetric, detT = Pf(T) 2 . But, as we see easily from the definition (3.1), 
|Pf(T)| < D{T). By Theorem 3.1, |Pf(T c )| = D(T), where T c has the canonical flux distribution. 
By Theorem A.l, the gauge can be chosen so that T c is antisymmetric. ■ 
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